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Abstract. This paper considers multi-dimensional affine processes with continuous sample paths. 
By analyzing the Riccati system, which is associated with affine processes via the transform formula, 
we fully characterize the regions of exponents in which exponential moments of a given process do 
not explode at any time or explode at a given time. In these two cases, we also compute the long- 
term growth rate and the explosion rate for exponential moments. These results provide a handle 
to study implied volatility asymptotics in models where log-returns of stock prices are described by 
affine processes whose exponential moments do not have an explicit formula. 



1. Introduction 

Since the introduction of the Black-Scholes model, many models have been de veloped to captur e 
empirical feat u res of financia l asset prices. Among them, models proposed by ICox et al,l 19851 ] . 



Heston 



1993 ] 



Vasicek 



19771 ] . and many others have been widely used by market participants 



because of their analytical tractability in derivative pricing in addition to their ability to reflect 
observed market phenomena. Later, common features of these models were unified to introduce 
the notion of affine processes with the so called canoni cal state spa c e. Th e general treatment of 
affine processes w ith this state space was conducted by Duffie et al. 200o| and later extended by 



Duffle et al. 



20031 ]. Q uite recently, studies on affine processes have been extended to more general 



Cuchiero et al 



201 ll ] and references therein. 



state spaces; see, e.g., 

A defining feature of affine processes is the logarithm of their Fourier transform is a li near 
function of the state. The regularity of affine processes, proved in iKeller-Ressel et al.1 20111 ] for 
affine processes on the canonical state space, links the aforementioned linear function to solutions 
to a system of (generalized) Riccati differential equations. This connection contributes to the 
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analytical tractability of affine processes and enables us to express the values of deri vative contr acts, 
whose underlying is modeled by affine processes, via the Fourier inversion formula (jLed [2004al ] and 
references therein). Moreover, this connection bridges affine processes and the theory of dynamical 
systems. Distributional properties of affine processes can be charac terized by dynamical behaviors o f 
solutions to the assoc i ated Riccati system . We r efer t he reader tolFilipovic and Maverhoferl 20091 ]. 



Glasserman and Kim 



20101 ]. 



Keller-Ressel 



201ll | . and 



Cuchiero et al 



20111 ] for recent developments 



in this direction. 

In this paper, we investigate long-term and blow-up behaviors of exponential moments of affine 
processes with the canonical state space K^xR™. We treat general multivariate affine processes with 
continuous sample paths, so called affine diffusions. This restriction of af fine processes to diffusions 
is imp osed because its transform formula has been well understood in iFilipovic and Maverhofer 
2009| . Currently the tra nsform formula for affine processes with jumps is being studied; see 



Spreii and Veermanl 20101 ] . The generalization of our results to affine processes with jumps is 
left as future studies. By focusing on affine diffusions, we are able to find sharp answers to the 
following two questions: 

Ql: Given an affine diffusion X, what are all possible vectors u such that Eexp(it T Xt) < 00 for 
any T > 0? For such a vector u, does the long-term growth rate Huit^oo T~ l log{E exp(-u T Xr)} 
exist? 

Q2: For a fixed T > 0, what are all possible vectors u such that Eexp(u T X$) < 00 for any 
S < T? For a vector u such that Eexp(u T Xs) is finite for all S < T but infinity for S = T, 
does the blow-up rate lim.sfr(T — S) log{Eexp(-u T X5)} exist? 

These questions are motivated by practical applications explained in the next paragraph, but they 
are also mathematically interesting. By focusing on a class of affine diffusions with some hierarchical 
structure between components (see Assumption [1]) , we provide complete answers to Ql and Q2. (It 
should be noted that this class contains virtually all affine diffusions with the canonical state space in 
financial modeling.) The set of u such that Eexp(w T Xr) < 00 for any T > is characterized via the 
disjoint union of stable sets for equilibrium points of the Riccati system in Theorem 13. 4[ Moreover, 
the growth rate of exponential moments is identified in Corollary 13.61 Working with a transformed 



Riccati system, similar answers to Q2 are proyided in T 
results are extensions of 



Glasserman and Kim 



10] 



reoreni 13.91 an d Co rollary 13.101 These 



2010 and 



Keller-Ressell 20111 ] to affine diffusions 



with arbitrary dimension. These findings not only help numerically identify sets of vectors in Ql 
and Q2, they also characterize large-time asymptotics and explosion phenomenon of exponential 
moments of multi-dimensional affine diffusions. 

For the past several years, large-time asymptotics and explosion phenomena of stock price mo- 
ments have attracted considerable attention because of their cl ose connecti on to implied volatility 
asymptotics. By approximating long-term stock price moments, 



Lewis 



20001 ] derived an asymptotic 
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formula for the i mplied volatility at large maturities in the fixed-strike regime under the Heston 
model. Recently, iForde and Jacquier |20lJ ] obtained similar implied volatility asymptotics for the 
Heston model in a regime where the log-moneyness is proportional to the maturity. The first step 
in their analysis is t o stu dy the long-term behaviors of stock price moments (see Theorem 2.1 in 



Forde and Jacquierl 201 11 ] ) . On the other hand, it is well known that the explosion of certain mo- 



ments of stock prices a t fixed time T is related to the im p lied v olatilities at extreme strikes with 
option maturity T; see 



Lcc 



2004bl ] and 



Benaim and Friz 



20081 ] for extensions. For example, an 



upper bound on the asymptotic slope of implied volatilities of deep-out-of-money options is found to 
be a function of the critical exponent p* = sup{p | ES^, +1 < 00}. Such asymptotic values of implied 
volatilities are informational in extrapolating smile curves and in calibrating u nderlying models 
to ma rket prices. More details abou t this practical usage can be found in, e.g., iBenaim and Friz 
20081 ] and lForde and Jacquierl [20111 ]. When the stock price log-return is modeled by an affine diffu- 
sion, results in this paper help to identify implied volatility asymptotics for large-time-to-maturity, 
deep-out-of-money or deep-in-the-money options; see Section 13.31 and three examples in Section |H 

The paper is structured as follows. In Section [2j we review basic concepts of affine diffusions 
and their canonical representations. We present our main results in Sectional Then, three multi- 
dimensional examples are presented to illustrate our findings in Section [4[ Analysis on the Riccati 
system and proofs of main results are developed in Sections [5] and (6[ Finally, Section [7] concludes. 

Before we move on, let us introduce some notational conventions which will be used throughout 
the paper. 

• For a vector x in a Euclidean space, |x| means its Euclidean norm regardless of dimension. 

• If x, y are of the same dimension then x < y if and only if Xi < yi for each component. And 
x ■ y represents the Euclidean inner product between x and y. 

• For a vector in M m+n or a matrix in ]j( m + n ) x ( m + n ) ) we denote the first m entries of the 
vector or m x m entries of the matrix by the superscript V, and the last n entries of the 
vector ornxn entries of the matrix by the superscript T>. 

• By , where x £ M m and X C {1, . . . , m}, we mean a vector of which i-th entry is 

• For matrices, diag(x) for x 6 M m is the m x m diagonal matrix with (x±, . . . ,x m ) as its 
diagonal entries, and diagx(x) with I C {1, . . . , m} is the mx m diagonal matrix such that 
its i-th diagonal entry is Xjljgx- Ik 1S t ne k x k identity matrix. 

• For a set A in Euclidean space, A° is its interior and A c is its complement. 



2. Affine Diffusions on Canonical State Space 

Let us recall affine diffusions and their canonical representation in this section. Given b : x 
W 2 — > M. d and a : M.™ x M. n — > M. dxd for some nonnegative integers m, n and d = m + n, we consider 
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the following stochastic differential equation (SDE) on a probability space (fi, (J~t)teR+i ^) : 

dY t = b(Y t ) dt + a{Y t ) dW t , Y = y, 

where W is a d-dimensional standard Brownian motion and y £ RT x R n . The above SDE admits 
a unique solution when b and a are of affine type and s atisfy admissible constraints introduced 
below (see Theorem 8.1 in lFilipovic and Mayerhoferi [20091 ] ). The state space x W 1 is called the 
canonical state space. In financial applications, the first m components of Y, which usually model 
volatility processes, are called volatility state variables; while the other n components of Y are 
called dependent state variables. In this case, Y v models the volatility variables and Y v describes 
the dependent variables. 

We say that Y is an affine process if there exist C- and devalued functions <f> and tp such that 



(2.1) 



E 



exp [u Y T ) \F t 



exp (V(T -t,u) + V(T - t, u) T Y t ) , 

for all u £ iM. d , t < T, and y £ R™ x R n . This specification implies that t he diffusion matrix a(y) := 
o~(v)& (y) T and the drift b(y) are both affine functions (see Theorem 2.2 in lFilipovic and Mayer hofei 



20091 ] ). i.e., 



a(y) = a + Vi a ^ b (v) =b+ ^2 ='■ b + B Vi 



i=l 



i=i 



for some a, on £ R dxd and column vectors b, /?,; £ R d with -B := (/3i ■ ■ ■ £ R dxd . Moreover, regu- 
larity of affine processes proved by iKeller-Ressel et all [20111 ] ensures that </> and tp = (ipi, ■ ■ ■ ,ipj) 



satisfy the following system of Riccati differential equations: 



d t (j>(t,u) 
dtipi(t,u) 



4>(t,u) T ai)(t,u) + b T ij(t,u), <f>(Q,u) = 0, 

^(t,u) J aitp(t,u) + (3jtp(t,u), ip(0,u) = u, for 1 < i < d. 



To ensure that Y is an affine process on the state space M.™ x R n , we impose the following 
admis sible constraints on parameters a,ai,b, and (see Theorem 3.2 in iFilipovic and Maverhofei 



20091 ]'): 



(i) a, o.i are symmetric positive semi-definite, and a m +i 

' o o \ / aA 



(ii) a 



a 



— C^m+n — 0, 

where Cj £ R, 5,, £ M mxm is the zero matrix except 



1 for the (i,i)-th entry, and w% £ R mxn has zero entries except the i-th row, 
( B v \ 

(iii) b £ WV; x R n , i? = ^ , and i? v has nonnegative off-diagonal elements. 

Under these constraints, the transformation formula (|2.ip is extendable to real dimensions. 



Theorem 2.1 (jFilipovic and Maverhoferl 20091 ]). Suppose that Y is an affine process with admis- 
sible parameters. Then, the transform formula (|2.1j) holds true for u £ R^ as long as either side of 
the formula is finite. 
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In this paper, we focus on the following class of affine diffusions: 
Assumption 1. B v is triangular (say, upper triangular) with strictly negative eigenvalues. 

The upper triangular shape of B v imposes a hierarchical dependence structure between all volatil- 
ity state variables. This hierarchical structure is commonly assumed in many financial models (see 
Section U] for several examples). In these models, different volatility state variables are usually used 
to model volatility processes on different time scales. Strictly negative eigenvalues imply that Y v 
is mean-reverting, which is a natural property of volatility processes. 

To facilitate our an alysis on this class of affine diffus ions, we consider their canonical representa- 
tions (see Section 7 in lFilipovic and Maverhoferl 20091 ] ). Given a linear transform A : M.™ x M. n — > 
R™ x R n , the process X := AY has the following dynamics: 

dX t = {b + BX t ) dt + a{X t ) dW t , X = AY , 

where b = Ab, B = ABA~ 1 , and <j(x) = Aa(A x). The transformed diffusion matrix is 

d m+n 

a(x) = a(x)a(x) T = AaA T + ^(A~ 1 x) i AQ i A T =: a + ^ Xi&f. 

i=l 1 

Actually, one can find a special A S M. dxd with diagonal A v , such that the diffusion matrix of X 
has the following canonical form 

diag x (x) \ 

v 7T + Y1T x i 7T i J ' 

where X is a subset of {1, • • • , m}, and 7Tj, < i < m, are some symmetric positive semi- 
definite matrices in M nxn Mo reover, the parameters of X are admissible. (See Lemma 7.1 in 



d(x) 



Filipovic and Maverhoferl [20091].) Note that B v is still upper triangular since A v is diagonal. 



Moreover, to exclude trivial cases where Y v have deterministic dynamics, we assume that T is 
non-empty. 

In the canonical version, the Riccati system reads 

d t (f>(t,u) = ^ v (t,u) T Tr ^ V (t,u)+b T i;(t,u), 

dti/>i(t,u) = ^M^ufhez + ^it^firi^i^^+pJ^u), l<i<m, 
dtifji(t,u) = fa ip(t,u), m + l<i<m + n, 

with initial conditions <p(0,u) = and tp(0,u) = u. Note that the first equation is easy to solve 
once we know ip, hence we focus on equations for ip and write them succinctly as follows: 

y = f(y,z), y(o) = v, 

z = A v z, z{0) = w. 
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Here u = (v, w) with v eM. m and w G R n , / = (/i, • • • , f m ) T with 



9{z) 



fi(y,z) :=\yfl 


\ 


' Z T 7Ti Z ' 


1 




"~ 2 


^ z T vr m z y 


/ 



fc=l 



+ ^ c z and 





Assumption Q] implies that is a lower triangular matrix. Hence (yi,--- ,yi) in (|2.2p is an 
autonomous system for each i G {1, • • • ,m} when w G Ker Moreover, has strictly negative 
eigenvalues with nonnegative off-diagonal elements, whence —A v is a nonsingular M-matrix (see 
Definition IA.3|) . Now the transform formula reads 



(2.3) 



E 



exp(u T X T )\T t ] = exp (l(T - t) + y(T - t) • X f v + z(T - t) • ) 



where /(•) := (1/2) / z(s) T vr 2;(s) ds + f Q b v ■ y(s) ds + J Q h v ■ z(s) ds and (y, z) solves (p^D . 

In financial applications, the discounted stock price, say S, is usually modeled by an affine process 
X via S. = exp(^ T X.) for some 6 G M. d . Then, S being a martingale under (a risk neutral measure) 



P implies that 9 V G Ker^l^ (see (|3.5p ). Therefore, in this paper, we always choose the initial 
condition for the second equation in (I2.2p to be z(0) = w G Ker A v . Hence, z(t) = w for any t > 0, 
and the first equation in (|2.2p reads 



(Ric-V) y = f(y,w), y(0)=v. 

We call v G M m an equilibrium point of (|Ric-Vp if /(v, io) = 0. 

Remark 2.2. It is of potential mathematical interest to consider models without Assumption [TJ 
However, in such cases, even identifying all equilibrium points of (|Ric-Vp becomes a nontrivial task 
as we need to solve a system of coupled algebraic equations. Still, there is one case where some of 
the results in this paper can be obtained to some extent, and this is when A v is invertible. We 
refer the reader to 



Kim 



20101 ] for details. 



3. Main Results 

In this section, we present our main results whose proofs are deferred to Sections [5] and [6l In 
Section EU we look for all u G M m x Ker A v such that E[exp(n T A T )] is finite for all T > 0. In 
Section [3.21 we characterize all u G W 71 x Ker A v such that E[exp(ti T Xg)] is finite for all S before 
a given T. Applications of these characterizations to financial modelings are given in Section 13.31 



ON EXPONENTIAL MOMENTS OF AFFINE DIFFUSIONS 7 

3.1. Long-term behaviors. Our first result identifies u G M. m x Ker A v such that E[exp(u T XT)] 
is finite for all T > 0. Thanks to Theorem 12. 1^ the problem is equivalent to finding every initial 
condition v G M m such that the solution y to (|Ric-Vp does not blow up in finite time. To this 
end, let us first classify equilibrium points of (|Ric-Vp into several different types, each of which 



tells us abou t qua 



itat ive behavior s of solutions in a neighborhood of an equilibrium point. See 



Chiang et al.l 19881 ] or iPerkol 20011 ] for more backgrounds. 



Definition 3.1. An equilibrium point v G K m of (|Ric-Vp is stable if for each e > there exists 
5 > such that \\y(t) — u\\ < e for all t > whenever ||y(0) — v\\ < 5. It is asymptotically stable if 
it is stable and lim^oo y(t) = v. Otherwise, v is unstable. Also, if all eigenvalues of the Jacobian 
Df(v) = j4 v + diagxiy) of / at v have nonzero real parts, then v is hyperbolic. 

The following result identifies all asymptotically stable equilibrium points for (|Ric-Vp . 

Lemma 3.2. There exists a nonempty closed convex set 2) C Ker A v with the following properties. 
First of all, for each w G 2), there are at most 2^ equilibrium points for (IRic-Vp . Second, for each 
w G D°, T](w) = (r]i(w), • • • , r] m (w)), where 



-An - Aa\ - 2 MJi AikVkH + 9i(w)) if i G 1 

(3.1) rji{w) := { V _ / \ 

-Ail 1 {£l=i A ikm(w) +gi(w)j if i G {1, • • • ,m}\l 

is hyperbolic and it is the unique asymptotically stable equilibrium point. All other equilibrium points 
are unstable, while at least one of them is hyperbolic. Lastly, there is no equilibrium point when 
w G KerA v flS c . 

The construction of D is explicit (see (|5.ip below). Moreover, rj(w) can be determined sequentially 
from i = 1 to i = m since A v is lower triangular with strictly negative diagonal entries. Now in 
order to connect the long-term behavior of solution trajectories to equilibrium points, we introduce 
the following notion. 

Definition 3.3. Given an equilibrium point v of (|Ric-Vp . its stable set is 

W*(w) := \ v G M' m | limy(i) = v where y(t) solves (|Ric-VD I . 

When v = rj(w), we write W£(w) as S(w) and call it the stable region of (jRic- Vp . 

Another related object is the set of initial conditions for (|2.2p such that its solution trajectory 
does not explode in finite time: 

Soo : = { n = ( v ,w) £i m x Ker A v \ \y(t)\ < oo for all t G K+, where y solves (jRic-Vp } . 

For each w, S QO (w) is the section of 5oo, i.e. S ao (w) := {v G M m | (v, w) G 5oo}. We are now ready to 
state our first main result, which provides a decomposition of Soq. The interior of 5^ is the disjoint 
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union of stable regions S(w) for all w G 3°, the boundary of Sqq consists of two components: 1. 
disjoint union of all stable sets of nonstable equilibria v for each w G 2. disjoint union of 
5oo(w) for each w G dD. In all of our statements, the topology is the relative Euclidean topology 
of R m x Ker^4 D . 

Theorem 3.4. The interior and the boundary of 5^ have the following decompositions in R m x 
Ker A v : 

i) = U w€ ^S(w) x {w}. 

ii) dSoo = flJ«,6£i°U^(«,)^H x U (u„,eas) 5 oo(w) x {w}^J , where v is chosen 
from equilibrium points of (|Ric-V|) . Moreover, for each w G dD, there exists a nonempty 
set M C {1, • • • , m} suc/i that the set {vm \ v G e>oo(i(;)} is i/ie stable set of the following 
system 

Vi = -jVi^iax + -4ifcyfe + fori€M, 
keM 

which admits a unique equilibrium point. Here v_m := (v^, . . . ,Vi k ) if M = . . . ,ik} C 
{1, . . . ,m}. 

In some special cases, the description of dS^ becomes succinct. A hyperbolic equilibrium point 
v of (|Ric-Vp is of type k if it admits k eigenvalues with positive real parts in its Jacobian matrix. 
It is a standard result in dynamical systems theory that the stable set of an equilibrium point 
v G W n of type A; is a smooth manifold of dimension m — k. If the system of interest has hyperbolic 
equilibrium points only, then the description of the (m — l)-dimensional object dSoo(w) for w G 3° 
does not need the stable sets for equilibrium points of type k > 1, because these stable sets have 
(m — l)-dimensional Lesbegue measure zero. 

Corollary 3.5. Suppose that A v is invertible and every equilibrium point of (|Ric-Vp is hyperbolic. 
Then, dSoo is given by [J u W„(0) x {0} except a set of (m — 1)- dimensional Lesbegue measure zero. 
Here, v is chosen from hyperbolic equilibria of type 1 and W*(0) is a smooth manifold of dimension 
m — 1. 



Going back to the afhne diffusion X, the characterization of 5oo, together with Theorem 12.1 
helps to identify the long run behavior of its exponential moments. 

Corollary 3.6. The following statements are equivalent: 

i) E [exp (u T X T )] is finite for all T > and Xq G M.™ x M. n . 

ii) u G 5oo- 

Moreover, when either of these statements holds true, 

(in\ v 1 ^ ml ( Tv \1 l (-D\T T>,rD D.J b v ■ i]^), if u G 5^; 
(3.2) lrni -logE exp(u X T ) = -(u ) ir u + b ■ u + < t 
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where v is some unstable equilibrium point of (|Ric-Vp . 



9 



These findings c onnect to existing results in two ways. First, it g eneralizes cha r acter izations in 



Glasserman and Kim 



20101 ] and Theorem 3.4 in 



Keller-Ressel 



20111 ] to multi- 



Proposition 5.2 of 

dimensional affine diffusions. In these two papers, similar charact erizations on exponential mo 
ments are obtained in the canonical affine term structure model of 



Dai and Singleton 



30] 



2000] and 



2-dimensional affine sto chastic volati 
ments in Theorem 3.4 of Keller-Ressel 



ity m odels, respectively. Second, following the same argu- 



201 ll ] . Corollary 13.61 shows a certain similarity between large 



time moment generating functions of X and a Levy process whose characteristic exponent is given 
by the right hand side of (|3.2p . 

3.2. Blow-up behaviors. Given T > 0, our second result identifies u G M. m x Ker A v such that 
E [exp (u T X5)] is finite for any S < T. To this end, let us first define the blow-up time for solutions 
to (IRic-Vp . 

Definition 3.7. For the initial condition u E W 71 x KerA^, the blow-up time T*{u) of a solution 
y to (|Ric-Vp is the first time t* such that \mi t _+ t * \y{t)\ = oo. 

The following result, whose proof is deferred to Section El ensures the continuity of u i— > T*(u). 

Lemma 3.8. The blow-up time T*(-) is continuous on the set V := {u \ T*(u) < oo}. 

Similar to 5oo in the last subsection, we define the set of initial conditions such that solutions to 
(|2.2p do not blow up before T: 

S T := {u = (v,w) El ffi x Ker A° \ \y(s)\ < oo,Vs < T, where y solves (|Ric-Vp } . 

We also define St(w) as a section of St for fixed w £ Ker^". It is apparent that St{w) = 
{v | T*{u) > T where u = (v,w)} and that {St(w)} t>0 is a decreasing sequence of sets, yielding 
C\t>o^t(w) = Ht>o i v I T*(u) > T} = Soo(w). This observation and Lemma 13.81 combined indi- 
cates that both Soo(w) and St("w) are closed sets in R m . Moreover, iFilipovic and Maverhoferl [20091 ] 
showed that St(w) is a convex neighborhood of the origin in M. m . Similar conclusions hold for St 
and Soo in M m x Ker A v as well. On the other hand, it is not difficult to see from the definition of 
S T that S T = {u \ T*(u) > T}, hence dS T = {u \ T*{u) = T}. 

In what follows, we will characterize St and its boundary via the stability analysis of a trans- 
formed version of (|Ric-Vp . Before we proceed, observe that St = {u \ E [exp (^X^)] < oo, \/S < T} 
from Theorem l2.11 Hence, the study of St and its boundary is equivalent to investigating the blow- 
up behaviors of exponential moments of X. 



Let us consider the following change of variables, inspired by iGorielvl [20011 ]: 

Xi (s) := e~ s yi (T(l - e~ s )) , i = l,...,rn, x m+1 (s) := e~ 
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Observe that if Xj blows up at some s* > 0, then yi blows up at T(l — e~ s ) < T. Therefore if 
yi explodes at T, then x% does not explode in finite time. In addition, if m explodes after T, then 
lim^oo x(s) = 0. Given y(0) = v G M. m , one checks that x satisfies the system of ODEs: 

rj-i 1 

(3.3) ±j = —x?I ie z - Xj + T A ik x k x m+ i + Tx 2 m+1 g(w), i = l,...,m, 

k=l 

with x m+ i = —x m+ i and the initial condition x(0) = (v, 1). We introduced the auxiliary component 
x m +i to ensure the system (|3.3p is autonomous. Then, we observe that the equilibrium points of 
(j3.3[) are given by v 1 with v[ = or 2/T, for i G I, and zero for all other indices. Also, every 
equilibrium point is hyperbolic, since the Jacobian at each equilibrium point has eigenvalues 1 or 
— 1. Furthermore, the origin is the unique asymptotically stable equilibrium point of the system. 
For each equilibrium point v' for (|3.3p . let us denote its first m components by v and define the 
following stable set for u: 

W°(w,T) := \ v G IR m | \imx(t) = v' where xit) solves d3T3j) \ . 

We are now ready to state our second main result, which characterizes the interior and the boundary 
of St as the disjoint unions of stable sets of equilibrium points for (13.3|) . 



Theorem 3.9. For each T > 0, the interior and the boundary of St have the following decompo- 
sitions in R m x Ker A v : 

S T = (J W$(w, T) x {w} and 8S T = [j j |J W°(w, T) x {w} j , 

w&Kcv A w&KcrA D / 

where v is chosen from the first m components of equilibrium points of (13.3|) . 

In the same spirit of Corollary 13. 6[ the blow-up behavior of exponential moments is identified as 
follows. 

Corollary 3.10. For each T > 0, the following statements are equivalent: 

i) E [exp {u T X s )] is finite for all S < T and Xq G x W 1 . 

ii) u £ St- 

If u G 8St, then 

(3.4) lim(T - S) logE \exp(u J 'X s j\ = Tu ■ XX , 

5tT L 

where v is the first m components of some unstable equilibrium point of (|3.3h . 
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3.3. Financial applications. Affine processes have been widely used to model the stock price dy- 
namics because of their analytical tractability in derivative pricing. In many models, the discounted 
stock price is represented by S. = ex.p(9 T X.) for some 9 S R rf . Let us assume that S is a martingale 
under a risk neutral measure P. Then, this assumption is equivalent to the following conditions: 

(3.5) is an equilibrium point of (H2J) and 1/2 (9 v ) T Tr 9 v + b T 9 = 0. 

In particular, 9 V £ KerA®. To prove (|3.5|) . we have from (|2.3|) that 



, E[St\ H 

1 = = exp 

Of. 



I(T-t) + (y(T-t)-9 v )-X? + (z(T-t)-9 v )-X 



r'D 



Vt < T and X u 



if and only if J(s) = and (y(s), z(s)) = (9 V , 9 V ) for any s E R+. Hence (|3.5|) is confirmed. 
Now, the long run behavior of stock prices in this model follows from Corollary 13.61 directly. 



Proposition 3.11. For A € R, the following statements are equivalent: 

i) E[S$] is finite for any T > and Xq £ R™ x R n . 

ii) X9 £ Soo. 

When either of the above statements holds true, the asymptotical growth rate of the stock price 
moment is given by 

(3.6) lim -logE[5 T ] = -A (9 ) tt 9 + \b ■ 9 + ■> , 

T^oo I l y b v ■ v, if X9 G oSoo 

where v is some unstable equilibrium point of (IRic-Vjl . 

The characterization above has implications on prices of securities with super-linear payoffs. 
Andersen and Piterbarel 20071 ] discuss possible unbounded prices of securities under two-factor 
affine or non-affine stochastic volatility models. Moreover, (|3.6p can be used for the large-time- 
to-maturity implied volatilities for European options in multi-dimensional affine models. These 
asymptotic form ulae fac i litate calib rating models to implied v olatility surfaces in practice and have 
been obtained in Lewie |2000 ] and Forde and Jacquier 2011 ] for one volatility factor models. We 



provide several examples of multi-dimensional volatility factor models in Section UJ 



Remark 3.12. In 



Forde and Jacquierl [20111 ] and 



Keller-Ressel 



20111 ] . a parametric constraint was 



imposed when the long-term growth rate was calculated. As we will see in Section HI this parameter 
constraint is equivalent to 9 being a stable equilibrium point. But, (|3.6p still holds even when 9 is 
unstable. 



The characterization of blow-up regions in Theorem 13.91 ties closely to the implied volatility 
asymptotics at extreme strikes for European options with fixed maturities. Let us denote a 2 (x,T) 
the implied volatility for a European option with strike K, maturity T and the log-moneyness 
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x = log(K/S ). 
(3.7) 



Lee 



2004bl ] proved that 

<? 2 {x,T) 
\x\/T 



lim sup 

x— >oo 



c(p ), hmsup 

X — > — oo l-^l/ 



where p* = sup{p > | E[Sf, + ] < oo}, q* = sup{g > | E[5 T 9 ] < oo}, and c(cc) = 2 -4(\/x 2 + x-x). 



lere p* and q* are called critical exponents. This result was extended later by 



Benaim and Friz 



20081 ] . where the limit superiors in (13. 7p are replaced by limits. These asymptotic values of implied 



volatilities at extreme strikes have been found to be useful for extrapolation of smile curves (see 



Benaim and Friz 



20081 ] ). It is then vital to calculate critical exponents of underlying models in 



order to apply aforementioned connections to implied volatility asymptotics. 

In the model where the logarithm of the discounted stock price is 9 T X, critical exponents can 
be identified by looking at dSx- Before we proceed, let us first re-define critical exponents because 
they depend on the initial condition Xq in our multi-dimensional setting. For example, consider a 
case where each component of Xq is zero as long as the corresponding component of y blows up at 
T. Then, the blow-up time of exponential moments would not be equal to that of y. Thus, we set 
p* as follows: 



p* := sup{p|E[^ +1 ] < oo,VX } 



sup |p I E 



exp((p + l)9 T X T )j <oo,VX } 



= sup{p| \y(T)\ < oo with (y(0), z(0)) = {p + 1)9} = sup {p | T*((p + 1)6) > T} , 

where the third equality follows from fj2.3[) and Theorem 12. 1[ the fourth equality holds since T*(-9) 
is nonincreasing (see Lemma 16.21 below) . We also note that p* > because the martingale property 



of S implies that E[exp(6> T X T )] = exp(9 Xq) < oo for any T. Similarly, we re-define q* := 
SU P {q I T*{— qO) > T}. Now it follows from the continuity of T* (see Lemma I3.8|) that T*((p* + 
1)9) = T*(—q*9) = T. Hence, the intersections of a line passing through the origin and 9 with dSx 
yield critical exponents. 



4. Examples 

4.1. Heston model. Let us start with the Heston model which is a prominent example of two- 
dimensional affine stochastic volatility model. The dynamics is determined by the SDE: 

«-((t)-(^:H* + ^(;^H 

where Wj is a standard two-dimensional Brownian motion. The discounted stock price is modeled 
by St = exp(Y^' D ), the variance process is described by Y v , and the diffusion matrix is a{y) = 

yi ( ^ | . Choose A = ( ^ j . The canonical version X = AY has dynamics dXt = 

V ap 1 I \ -p/a 1 J 
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Figure 4.1. Decomposition of Soo for the Heston model with k = 1, o = 0.4, and 
p = -0.2. 

+ BX t )dt + o(X t )dW t with 




B 



-K 

npo -a 2 / 2 




p 



yi, 



and the initial condition is given by X$ = (Vq/o 2 , —pVo/o + log So) where So, Vq are the initial 
stock and variance levels. Moreover, log(SV) = (p&, 1) ■ Xt, thus 6 = (pa, 1) T . 

For any initial condition (y(0), z(0)) = (v,w) of (|2.2p . z(t) = w and solves y = (l/2)y 2 — 
Ky + g(w) where g(w) = a 2 (l — p 2 )w 2 /2 + (npa — a 2 /2)w. It is clear that Ker^" = R. In addition, 
2) = \w | 2g[w) < k 2 } and 3}° = {u; | 2g[w) < k 2 }, following the definitions in Section [5.11 For each 
w € 3°, the equation for y admits two equilibrium points: L(w) = K — \/ k 2 — 2g(w) and U(w) = 
k+ y k 2 — 2g(w). The former is asymptotically stable and hyperbolic, and the latter is unstable and 
hyperbolic as the Jacobian Df(w) = —^/n 2 — 2g{w) < (\/ k 2 — 2g(w) > 0) at L(w) (U(w) resp.). 
Therefore, the stable region is S(w) = (— oo, U(w)) hence <S^ = {(v,w) \ w £ 55°, v < U(w)}. The 
boundary dS^ is readily obtained as well. These sets are illustrated in the left panel in Figure 14.11 
for one set of parameters. 



Let us now comment on the parametric constraint k > op in iForde and Jacquierl 20111 ] . This 
constraint is actually the necessary and sufficient condition for 9 being a stable equilibrium point, 
in particular, 1 G 2)°. Indeed, recall that 6 = (op, 1) T is an equilibrium point, hence 9 G 6*oo and in 
particular 1 G D. When k = op, k 2 — 2g(l) = 0, then 1 G <9£>. When k < op, U(l) = k+\k — op\ = 
op hence 9 G dSoo and unstable. 

Under the parametric constraint in the last paragraph, extending 9 in both directions until it 
reaches dS^, we obtain p±9 with p + > 1 and p_ < 0. This is illustrated in the right panel of 
Figure I4~T1 Now it follows from Theorem 13.41 that E[5j.] is finite for any T > and p G [p_,p+]. 
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FIGURE 4.2. 8St for the Heston model with k = 1, a = 0.4, and p = -0.2. 



Actually, Proposition 13.111 implies that 

A(p) := lim — logEexp (p9 T Xt 



KLpp Kip 

-p+^rL(p), p e (p-,p+). 



This result coincides with Theorem 2.1 in 



a 



Forde and Jacquierl 201 If ] where the authors continue to 



prove the essential smoothness of A(-) and derive for mulae for the large-time-to-maturity implied 
volatilities. On the other hand, in iKeller-Ressell 20111 ] . the author provides the same formula under 



the same constraint and argues that the price process gets close to a NIG Levy model, in terms of 
marginal distributions. 

Figure 14.21 shows 8St for several different T values. From this figure, we can identify critical 
exponents p* and q* for fixed T as the first positive numbers p and q such that {p + 1)0 and — q9 
belong to 8St- We can also clearly see the convergence of St to 5oo as T — > oo. From the viewpoint 
of Theorem 13.91 u £ 9St translates into the condition that (v, 1) is the initial condition x(0) such 
that lim^oo x(s) = (2/T, 0) where x(-) solves (|3.3p . Therefore, as in Figure [4T31 one finds v such that 
(v, 1) is on the boundary of the stable set of (|3.3[) . yielding (v, 1) S 8St- However, we also note that 
the Heston model admits a closed form formula for 8St by which implied volatilities at extreme 
strikes can be calculated. This theorem, nevertheless, provides one method of accomplishing the 
same task even when such a closed form formula is not available. 



4.2. Double st ochastic volati lity model. The following 3-dimensional stochastic volatility model 
was proposed in iGatherall [20081 ]: 



dV t = Kxfyi -V t )dt + V t a dZ}, 



dV t 
dS t 



K 2 {v-V t )dt + V^dZl 



StJVtdZ, 



t ■ 



where Z% is a correlated 3-dimensional Brownian motion and k±, <p are strictly positive constants. 
In this model, V models the high-frequency variance and V' represents the low-frequency variance. 
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Figure 4.3. The stable boundary of (f3T5|l with w = 5 and T = 3. 



For this model to be an affine diffusion, it is necessarily that a = /3 = 1/2 and Z 2 is independent of 
Z 1 and Z 3 . Denote the correlation between Z 1 and Z s as p G (— 1, 1). We assume that ki > «2 so 
that the mean reverting speed of high-frequency variance is larger than its low-frequency analogue. 
It is an easy matter to check that Yt = (Vt, V/, log St) satisfies 





( 


( 









( 




«i 





dY t = 










+ 














{ 


\ 





J 




{ 


-1/2 










Ft + 



/ 07 \ 

y 7 ^ 
V ^/{l- P 2 )Yl ) 

where W is a 3-dimensional standard Brownian motion. Now, applying an appropriate linear 
transform X = AY, we obtain from straightforward calculations that X satisfies dXt = (6 + 
BX t )dt + 6{X t )dW t with 





( 


> 




6 = 










{ 


o J 





-«1 





, a(x) 



xi 
x 2 








\ 



V (l-p 2 )xi J 



Kl \ 
-K2 

p«i — 1/2 —pK\ / 

and the initial condition is given by Xq = (Vo, Vq, — pVo + log So). In addition, log St = Yj, 
(0 1)A~ 1 X T , which implies that 9 = 0,0, 1). 
The associated Riccati system ()2.2[) is 



(4.1) 



yi 



V\ - Km +gi(z), 



+ KlVl +52(2), 



and z = with the initial condition (yi(0), 2/2(0), z(0)) 
(p«i — 0.5)2 and 52(2) = —pK\Z. Clearly, Ker^- = R. 
Section 15.11 

£{w) = {(ui,u 2 ) I (iti - ki) 2 < k 2 - 2#i (w), (u 2 - K 2 f 



(ui,U2,w). Here, gi(z) = 0.5(1 - p 2 )z 2 + 
. Consider the following sets defined in 



< ni 



2/ciiti - 2^2 (w)} 



25 = {w G Ker^l «(«;) / 0} = G Ker^ D | /ef - 2px(u;) > 0, k\ - 2g 2 {w) > 2Km(w)} 
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Figure 4.4. Region 5^ and S T for T = 0.5, T = 1 with K\ = 2,k 2 = 1. 



with r]i(w) := ki — y/«f — 2g\{w). Then, Proposition 13.111 gives us 



A(p) := Jim ^logE[S£] 



K2<pm(p) forpeD , 772 (p) := k 2 - - 2# 2 (» - 2/ci7/i(p). 



Figure 4.3 shows the sets <95oo and for some specific set of parameters. Two graphs are 
obtained by numerically solving (|4.ip and (|3.3[) . Theorems 13.41 and 13.91 help to identify dSoo and 
35t numerically, because it suffices to find stable sets for unstable equilibria of their associated 
systems. To be more specific, for dSoo, we first locate all unstable equilibrium points of (|4.ip 
for fixed w. Then for each unstable v, we solve (I4.ip backward in time to retrieve its stable set. 
Finally, <95oo is obtained by patching all dSoo(w) together. The right panel of Figure H3] shows 
8St for T = 0.5 and 1. It is produced similarly by working with (I3.3p . Even though we only show 
part of dSx in the right panel, it is understood that the critical exponents p* + 1 and — q* are the 
intersections of the w-axis and dSx- 

To derive the large-time-to-maturity implied volatilities in this model, we need the following 
parameter restriction: 

Ki > p. 

As we have seen in the previous subsection, this restriction ensures that is a stable equilibrium 
point of (|4,ip . As a result, [0,1] C 3°. On the other hand, one can check that A is essentially 
smooth in D. Indeed, 



A'(p) = KlK 2 <P X 



-P + ffi(p)/\Ai ~2gi(p) 



-2 - 2 92(p) ~ 2«i77i(p) 

and both terms under the square roots converge to zero as w — > dT>. Consequently, the Gartner-Ellis 
Theorem applies, and thus {(log St — log So) /t] satisfies the Large Deviation Principle under P with 
the rate function (Legen dre transform) A*(x) = sup pG ^,{xp— A(p)} . It then follows from Proposition 
4.1.3 in I Jacquierl 20101 ] that the large-time-to-maturity implied volatility for the European option 
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with maturity T and strike price K(T) = Soexp(xT) is 
(4.2) 



a 2 (x, oo) := lim a 2 (x,T) = 2 [2A*(x) -x + 2 (l{ xe ( x * - l{ x <£{x*, *•*)}) A*(x) 2 ~ xA*(x) 



T->-oo 

where x* := A'(0) = — ip/2 and x* := A'(l) = 0.5ftiy?/(«4 — p). In particular, the large-time-to- 
maturity implied volatility for at-the-money European option is 

2/n a a K*fn\ «»/ n _ 1 - 2 P K X + \p\ y/l + Ak{ - A P K X 



<H(0,oo) = 8A*(0) = -8A(po), po 



2(1 -p2) 

Here, po is chosen so that A'(po) = 0- On the other hand, we can also obtain the leading order ex- 
pansion of cr(x, oo) when x is close to 0. This leading order expansion provides us information on the 
implied volatility asymptotics for fixed strikes. This is because one can choose x = T~ l log(K / Sq) 
for the fixed strike K. When T is large, x is close to zero. To obtain this expansion, one first 
observes that A*(x) = A*(0) + x(A*)'(0) + o(x) = A*(0) + p x + o(x) where the second identity 
follows from (A*)'(0) = (A')~ 1 (0) = po- Plugging the previous expansion into (14.2p . we get 

a 2 (x,oo) = 8A*(0) +2(2p - l)y/x + o(y/x). 
4.3. C ascading affine d iffusions. To take the multi-frequency aspect of interest rates into ac- 



count, |Calyet_et_aL| [20101 ] consider a model in which the interest rate has the dynamics that depends 
on several latent variables with high to low frequencies. The authors also suggest a multi-frequency 
stochastic volatility model for equity option pricing. In this subsection, we consider a specific form 
of such cascading volatility models which are also affine diffusions. 

In this model, Y v has no restriction other than the admissibility constraints on its parameters, 
while Y v follows 



dYl = k5 1 ~ (Y t l+1 - Y t l ) dt + aJYidWl i = 1, . . . ,m - 1 



dY t m = kS™- 1 (ip - Y t m ) dt + a^fYfdWr 

with < 5 < 1 and positive constants k, a, and (p. This model proposes that the volatility process 
Y 1 depends on many latent variables that have slower mean reversion speeds. Then, the process 

X = AY with A = | m | makes X a canonical version of Y. The associated Riccati 

V L, 



system (|2.2p reads 

Vi = \v 2 i + K {^~ 2 yi-A>i - $ l ~ 1 y l ) + ^z T ir iZ + (A c z) i , i = l,...,m, 

and i = A v ' z with the initial condition (y(0), z(0)) = «£ M. m x Ker A v . 

Let us consider a simple case where A v is invertible. Then, we easily see that Ker A v = {0 n } with 
n-dimensional zero vector n and that the equilibrium points are the origin, i.e., ry(O n ) = € M m , 
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and v = (0, . . . , 0, 2K5" 1 " 1 ). Moreover, the Jacobian matrices of equilibrium points are given by 
f J(0) M = -d^ 1 \ J( v ) iti = -<5 i - 1 I i<m + 5 m -% 



J(0) i+lti = 6*- 1 1 J{ri)i - ■ - v- 1 



Therefore, v is a hyperbolic equilibrium point and its stable submanifold W£ has dimension m — 1. 
Theorem 13. 4| then, implies that 5^ = <S(0 n ) x {0 n } and <95oo = W*(0 n ) x {0 n }. As a consequence, 
we obtain from Corollary 13.61 that 



1 r / \i //c<5 m_1 \ 2 

lim - logE exp ( u T X T ) = 2ip 

T-*oo T L V /J ycry 



for u G 95oo. As for the blow-up region, we have dSx = [y) u ^oW£(0,T)j x {0} where v is any 
m-dimensional vector whose entries are either or 2/T. 

Lastly, let us consider an equity model based on cascading affine diffusions. We just add one last 
process as follows: 

dY m+i = -iyj-dt + ^dW t m+1 

with a standard Brownian motion independent of (W 1 , . . . , W m ), and set St = exp(Y^ n+1 ) so that 
KSf = Kexp(p9 ■ Xt) where 6 = (0, ...,0,1). This can be understood as a general version of 
the Heston model but with no correlation between Brownian parts for the sake of simplicity of 
exposition. The associated Riccati system can be calculated accordingly. From Proposition 13.111 
we have A(p) := limy T" 1 log E5j. = Hnp5 m ~ l T) m (jp) / 'a 2 if pd G and u m if p9 G <95oo where v is 
some unstable equilibrium point. The stable equilibrium point of the above quadratic system can 
be found iteratively as follows: 



mi?) = K - V^ 2 - cr 2 p(p- 1), 



Vi (p) = kS^ 1 - ^kW-V -2K^- 2 r ?i „i(p), i = 2, . . . , m 

where p belongs to some interval, say [a,b], so that all the square root terms are well-defined. 
Such conditions read p G [p_,p+] with p± = (^a ± \/a 2 + 4k 2 ^ /(2a), and r)i-i(p) < k5 1 /2 for 
i = 2, . . . , m. It is useful to check that this interval gets strictly smaller as i increases, which we 
leave it as a simple exercise. Moreover, we get 



A'(p) 



a 



2 



X 



^ K 2 5 2{m-1) - 2K5 m - 2 7] m - 1 (p) 



At the boundary points of [a, b], the square root term converges to zero. Hence, |A'(p)| — > oo and 
thus A(p) is essentially smooth. By following the same arguments as in the previous subsection, we 
can obtain the implied volatility asymptotic formula at large-time-to-maturities. 



on exponential moments of affine diffusions 
5. Analysis of the long-term behavior 
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The long-term distributional properties of affine processes are determined by the long-term be- 
haviors of solutions for the associated Riccati system. Therefore, we shall first focus on equilibrium 
analysis of (|2.2p and prove Lemma 13.21 in Section [5. II and then characterize the long-term behavior 
of its solutions in Section 15.21 Finally, Theorem 13.41 and Corollary 13.61 are proven at the end of this 
section. 



5.1. S table equilibrium points. Let us start with some definitions inspired by iKeller-Ressel 



201 11 ] . Define the following two sets: 

<£ := ju G R d \ fi{u) < 0,V1 < i < mj and £° := ju G R d \ fi(u) < 0, V 1 < i < ?n| . 

They are sets of points on which all components of / are simultaneously (strictly) negative. It follows 
from the continuity of / that G; is closed (open), respectively. Moreover, they are convex thanks 
to the convexity of /. It is also clear that (£ is nonempty, since 6 f . Given w G KerA^, we 
define sections of £ and <£° as £(w) := {v £R m \ (v,w) G <£} and £°(w) := {v G R m \ (v,w) G <£°}. 
To identify all stable equilibrium points for (jRic-Vj) . we define 

(5.1) £> := {w G KerA v \ <£(w) / 0} and £° := {w G Ker A v \ £°(w) / 0} . 



It will be shown in Lemma 15.31 blow that £)° is indeed the interior of D. The first result below 
identifies the candidate stable equilibrium point for (|Ric-V|) . 

Lemma 5.1. Given w G T) , r](w), defined in (13. If) , is inside (£. Moreover, f(r](w),w) = and 
n(w) < v for any v G £(w). 

Proof. We utilize the lower triangular shape of A v and prove the statement by induction on i. This 
type of argument will be used repeatedly in our analysis. 

For i = 1, if i ^ X, f\ is linear with slope An < 0, then r)i(w) is chosen as the solution to 
/i(-, w) = 0. Clearly rj\{w) < v\ for any v G <B(w). If i G I, the quadratic equation fi(-,w) =0 has 
solution(s) because the graph of fi(-,w) has a nonempty intersection with M x R_. Then r]i(w) is 
chosen as the smaller of the two solutions (possibly the same) to the previous quadratic equation. 
It is also clear that r]i(w) < v\ for any v G £(w). 

Suppose now that the statement holds for k = 1, — 1. If i ^ X, it then follows from 
An < and A ik > for k / i that, for any v G £(w), rji(w) = -A^ 1 fefc=i AkVkiw) + gi(w)^j < 
—A~^ fX)fc=i Aik v k + 9i( w )j ^ v ii where the second inequality holds since fi(v,w) < 0. Now if 
i G X, notice that vf/2 + AuVi + YX=i A ikVk{w) + gi(w) < vf/2 + AuVi + X!=i Mk^k + 9i{w) < 0, 
for any v G <£(w). Then r]i(w), which is the smaller root of the quadratic function in Vi on the left 
side of above inequalities, must be less or equal to Vi. Hence the induction step is proved. Finally, 
r](w) G <B(w) is clear from the construction of rj{w). □ 
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Note that the Jacobian of / at r)(w) is A v + diagz{r]{w)) . A straightforward induction on the 
index i shows that w £ 3° if only only if A v + diagx(rj(w)) has strictly negative diagonals. In order 
to show that r](w) is indeed the unique asymptotically stable equilibrium point, let us present some 
topological properties of S3 and 3° in the next two lemmas. 

Lemma 5.2. The set S3 (T>°) is a nonempty closed (open) convex subset ofKerA v . Also, n is a 
convex function on 23 and i] £ C(23) n C 2 (S3°). 

Proof. Let us first show £ 23°. Indeed, for any v £ 2(0), f(v,0) = h v x + A v v < 0, and thus 
v > — (l/2)(A v )~ 1 Uj 2 ^ > 0, where the last two inequalities follow from the nonsingular M-matrix 
property of —A v (see Definition IA.3|) . Hence, Lemma 15 . 1 1 implies that n(0) must be 0. Notice that 
the Jacobian of /(-,0) at is A v , which has strictly negative diagonals. Therefore, £ 3° C 23. 

For the convexity of 23, it suffices to show that w := Xw + (1 — X)w £ 23 for any w, w £ 23 and 
A £ [0, 1]. Actually, since f(n(w),w) = f(r](w),w) = 0, it then follows from the convexity of / that 
f(Xt](w) + (1 — X)rj(w),w) < 0. This implies <£(w) ^ 0, hence w £ 23. Therefore we get n{w) from 
Lemma 15. II and n(w) < \rj(w) + (1 — \)r)(w) as rj(w) is the componentwise minimum of <B(w). 

The convexity of 23° is confirmed once we show w £ D° for any w, w £ D°. When w and w are in 
23°, both A v + diagx(i](w)) and A v + diagx(i](w)) have strictly negative diagonals. We deduce from 
the convexity of rj in D that A v + diagx(v(w)) < A (A v + diagx{i]{w)) + (1-A) (A v + diagx{rj(w))) . 
Strictly negative diagonal entries of the matrices on the right hand side imply that those of A v + 
diagx(i](w)) are also negative, which in turn implies w £ 2)°. 

That r\ is continuous is immediate from its construction in Lemma 15.11 Note that the square 
root term in rj(w) is nonzero for any w £ £>°, hence rj £ C 2 (D°) follows. The continuity of rj implies 
that D is closed. To see this, let us take a sequence {w n } C D such that lim n w n = w. Apparently, 
w £ Keryl 15 . Since rj is continuous, lirm^oo rf(w n ) = rj(w) which is defined in (13. ip . resulting in 
£(w) 7^ and thus w £ D. The openness of D° is obvious from the continuity of /. □ 

Lemma 5.3. If (£(io) ^ $, then <£°(\w) ^ for any A £ [0, 1). Therefore, D° is the interior of D. 

Proof. We begin with setting <£-i{w) := {v £ M m \ fi(v,w) < 0} and define (£° similarly. Certainly, 
£( w ) = C\iLi (similarly for <£?). 

Let us fix w £ S3. The proof is by an induction on i. When i = 1, if 1 ^ X, fi(-,w) is a linear 
function. Since An < 0, it is clear that (El(Xw) ^ for any A £ [0, 1). If 1 £ Z, fi(-,w) is quadratic, 
we denote the determinant of f±(-, Xw) as Ai(A) = — X 2 wTttiw — 2A(^4 c u;)i + A n , which is either 
a linear or quadratic function of A. It follows from <Bi(w) ^ that Ai(l) > 0. On the other 
hand, observe that Ai(0) = A\± > and w T 7Tiw > (tt\ is semi-positive definite), we then obtain 
Ai(A) > for any A £ [0, 1). This implies that <E\{Xw) / for any A £ [0, 1). 

Next, assuming that f|fc=i ^fc(Aw) / for any A £ [0, 1), we want to show that f|fc=i ^fc(Aw) / 
for any A £ [0, 1). For any k = 1, ■ • ■ , i — 1 and A £ [0, 1), rjj.(Xu)) is well defined by the induction 
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assumption and Lemma [5TT1 Consider v = (771 ( Aw;), •• • , rji-i(Xw), Wj, • • • ,v m ). Then, 



21 



1 i_1 1 

/»(«, Xw) = -vf I iez + A&Uj + r? fc (Au;) + -\ 2 w T iriW + A(A c w;)i 

fc=i 

is either a linear or quadratic function in Wj. As seen in the 7 = 1 case, the linear case is easy to 
handle. Hence we consider the quadratic case only. In this case, the determinant of f(-,Xw) is 
given by 

i-l 

Ai(A) = -A 2 w; T 7T;w; - 2A(A c w)i - 2 ^ ifc r/ fc (Au;) + 

fc=i 

Since Ajfe > 0, 7Tj is semi-positive definite, and 77 is convex (see Lemma |5.2|) . Aj(A) is concave in A. 
At A = 1, we know that 77(11;) exists and thus fi(v,w) = has a solution. Therefore, Aj(l) > 0. 
This observation together with the concavity of Aj(-) and Aj(0) = A\ > shows that Aj(A) > 
for any A £ [0, 1). Hence, Ha;=i &k(^ w ) ^ closing the induction. 

For the second statement, notice that for any w in the interior of T>, there exists r > 1 such that 
tw G 3). Then, the first statement implies Xtw S 3° for any A £ [0, 1). Hence, with A = 1/r, we 
see w £ 35°. Together with the openness of 23°, we conclude that 53° is the interior of D. □ 

Since A v is assumed to be lower triangular, the Jacobian Df(v) at a hyperbolic equilibrium point 
v has nonzero real eigenvalues. Moreover, it is well known that v is unstable if Df{y) has a positive 
eigenvalue, and v is asymptotically stable if all eigenvalues of Df(v) are negative. Consequently, a 
hyperbolic point v is (asymptotically) stable if and only if every eigenvalue of Df(i>) is negative. 
Now we are ready to prove Lemma [ 



Proof of Lemma \3.2l For w £ 33, all equilibrium points are constructed by solving f(v,w) = 
sequentially from index i = 1. Clearly, there are at most 2^ equilibrium points for (IRic-Vj) . 

When w £ 33°, we already observed that A v + diagx{i]{w)) has strictly negative diagonals. Thus, 
r](w) is hyperbolic and asymptotically stable. For the uniqueness, we simply note that, for another 
equilibrium point v, if we let i be the first index such that v\ > T]i(w), then by construction, it must 

be that i £ X and vi = — An + ^ Af { — 2 ^X^fc=\ ^ikVkiw) + gi{w)^ . Since the i-th diagonal entry 
of Df(v) is An + Vi > 0, hence v is unstable which contradicts with the choice of v. 

There exists at least one hyperbolic unstable equilibrium point. Take v such that its first I — 1 
components are equal to those of 77(77;) but vi = —An + Af t — 2 ^X)L=\ ^-ikVk( w ) + 9i(w)j where 
I = max {i \i £ X}. 

For w £ Ker A v n 33 c , if there exists an equilibrium point v for (|Ric-Vp . then (v, w) would be an 
equilibrium point for (|2.2|) . Therefore, v £ £(w) which contradicts to 77; £ 33 c . □ 
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Figure 5.1. Region <£°(w) and equilibrium points for a = 0.3,0.5,0.6. 



Example 5.4. For an intuitive understanding of the proof above, let us take a look at the following 
two dimensional system: 



1 



y 



f(y) = g 



yl 
y\ 



+ 



y, a > 0. 



-a 
1/2 -1 

In this example, I = {1,2} and m = 2. We set g(-) = for the simplicity of illustration. Therefore 
all sets blow are independent of w. Figure 15.11 shows the graphs of points that satisfy each of 
fi(y) = for three different values of a. The intersections of solid and dash lines are equilibrium 
points. The filled regions without boundaries correspond to their <£°(w) (for any w). We note that 
the number of equilibrium points changes as a varies. It is easy to check that the origin is hyperbolic 
whose Jacobian has negative eigenvalues as long as a > 0, or equivalently <B°(iu) ^ 0. At a = 0, 
<£°(w) becomes empty (also 2° = 0) and (0,0) is not hyperbolic. 

Observe that the point (1, 1) in the second panel is the only non-hyperbolic equilibrium point 
among all equilibria in three figures (and actually for all a > 0). Therefore, we can imagine that this 
kind of non-hyperbolic equilibrium case "rarely" happens. Inde ed, the hyperbolicity of equilibria is 
a generic property of dynamical systems (see, e.g., Section IV of lChiang et al.l |1988l ] and references 
therein) . 

□ 



5.2. Stable regions. After identifying the unique stable equilibrium point for (IRic- VP . we will 
study its associated stable region defined in Definition 13.31 When w £ £ > °, the stable region S(w) 
is an open neighborhood of r](w) diffeomorphic to W 71 (see IChiang et al.l 19881 ] ) . The next lemma 
gives a useful property on solution trajectories of (|Ric- Vp . 



Lemma 5.5. If the trajectory of (|Ric-Vp is bounded, then it converges to an equilibrium point. If 
the trajectory is unbounded, then it blows up in finite time. 



Proof. We use an induction starting from the index i = 1. But we present the induction step only, as 
i = 1 case is a specification of the general argument. Suppose that the trajectory of a solution y(t) 
is bounded in M m and that yk (t) for k = 1, . . . , i — 1 converges to the first i — 1 components of some 
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equilibrium point of (|Ric-Vp . If i ^ I, then yi(t) satisfies = e Aiit n/i(0) + J * e~ AiiS h(s)ds 
where h(s) = Y^k=i-^-ikUk{s) + 9i{w). Since An < and h(s) has a limit, say h(oo), it is easy to 
check that lim^y^t) = — h(oo)/An. 

If i E X, then we have a quadratic ODE = x(t) 2 /2 + /i(t) where x(i) = yi(t) + and 
/i(t) = Y^k=i AikUk(t) + gi(w) — A\j1. The induction assumption says that h(oo) exists. Since the 
trajectory of x(t) is bounded by assumption, Lemma |A. II implies that x(oo) := limtx(i) E M. m and 
x(oo) 2 /2 + h(oo) = 0. Then the limit of yi(t), together with (yi(i), . . . converges to the 

first i components of some equilibrium point for ([Ric-Vp . 

Let us move onto the second statement. Assume that i is the first index such that the trajectory 
of yi(t) is unbounded. It is necessarily that i E X. Let x be defined as above: x = x 2 /2 + h(t). By 
assumption, h{t) is bounded (actually, converges to h(oo) as we have seen in the last paragraph) and 
x(t) has unbounded trajectory. Therefore, we can find some time to and nonnegative constant C 
such that x(to) > y/2C and h(t) > — C. Let us consider a function x(t) which solves x = x 2 /2 — C for 
t > to with x(to) = x(to). Then, by the comparison theorem for scalar ODEs, we have x(t) > x(t) 
for all t > to- However, x(to) is greater than \J2C which is the unstable equilibrium point of the 
system for x{t). Elementary computations show that x(t) blows up in finite time, hence so does 
x{t). □ 

Remark 5.6. Following similar arguments as in the proof of Lemma l5.5( we can show that y{t) is 
bounded below for any initial conditions. 

Corollary 5.7. Let w € Ker A v n D c . Then, any solution y(t) to (IRic-Vj) blows up in finite time. 

Proof. If not, then the trajectory is bounded by the second statement of Lemma 15.51 Then the 
first statement of the previous lemma implies that it should converge to an equilibrium point, say 
v . This contradicts with the choice of w. □ 

We now take on the task of characterizing Soo(w) for w E D. Lemma 15.51 implies that S 0O (w) = 
\J u W£(w), where v ranges over all equilibria of (IRic-Vj) . Dealing with w E T>° and w E dD cases 
separately, we extract more information as presented in Lemmas 15.81 and | 



Lemma 5.8. Let w E D° . Then, Soo(w) = S(w) and <95oo(u>) = dS(w) = U^^(«,) W£(w) where 
is chosen from all equilibria of ()Ric-Vp . 



Proof. Part of the proof employs arguments in the proof of Theorem 4.1 of iKiml [20101 ] . but we 
include it here for completeness. Instead of (IRic-Vj) . it is more convenient to work with a slightly 
modified system on r(t) := y{t) — n(w): 

(5.2) f = ^r { z ) +Ar, A := A v + diagx(n(w)). 
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Since w G A has strictly negative diagonals. Additionally, A v is lower triangular with nonneg- 
ative off-diagonal entries, and therefore —A is a nonsingular M-matrix. It follows from Lemma 13.21 
that is the unique stable equilibrium point of (|5.2|) and its associated stable region S' is given by 

S' := \ r(0) | lim r(t) = where r solves (|5.2p > = S(w) — rj(w). 

Since S' is open, it contains a small neighborhood of 0. The set S'^ is defined analogously. Then 
it suffices to prove the corresponding statements for S' and S'^. 

To show the first statement, we begin with the observation that S 1 C S'^. Hence S' C S'^, since 
•S^j = 5oo(u>) — 7](w) is closed (see discussions after Lemma [378]) . In what follows, we will prove 
^oo \ 5' = by contradiction. Assume otherwise and pick r(0) € S'^ \ S'. It then follows from 
Lemma 15.51 that v := lim^oo r(t) is an equilibrium point. Since S' contains a neighborhood of 0, 
we can find 5 G (0,1) such that <5r(0) G dS'. Consider z(t) which solves (|5.2p with z(0) = 6r(Q). 
Recall that —A is a nonsingular M-matrix, and then we get that z(t) < 5r(t) for any t > from 
Lemma lA.4[ On the other hand, since z(0) G 55' C S'^, Lemma 15.51 implies that z(t) converges to 
an equilibrium, say u, of (15. 2p . Hence, v < bv. Moreover, these two equilibrium points v and v are 
nonzero; otherwise, r(t) or z(t) should have started from S' , which contradicts with their starting 
values outside S'. 

Note that v (also v) satisfies v = — (l/2)A~ 1 v^ > where the non-negativity follows from the 
property of a non-singular M-matrix —A. Moreover, < v < bv implies that 

2 x ~ 2 1 

from which we get —A(p — b 2 u) < 0. By multiplying both sides by — A -1 > 0, we obtain v < b 2 v. 
Repeating the same argument, we arrive at v < b 2k v for any integer k > 1. However b € (0, 1), 
hence v = 0. This contradicts with v ^ from the last paragraph. Therefore it is necessarily that 
5^ \ S 1 is empty. 

Since dS'^DS' = dS'nS' = 0, it is immediately seen that dS' C (Jj,_^ VF*' where the W,f' are the 
stable sets of each nonzero equilibrium point of (|5.2|) . To show the opposite inclusion, let us take 
an equilibrium point v ^ 0. Then, we claim that Wf)' C 95'. We prove this by contradiction. Let 
us assume that r(0) € W£'\dS' and lim^oo r(t) = z^. Then it is necessarily that r(0) ^ 5' because, 
otherwise, = lim^ r(t) = 0. Hence, r(0) ^ 5' and thus we can find <5 6 (0, 1) such that br(0) S 05'. 
By using the same argument as that in the last paragraph, we arrive at a contradiction. Hence 
dS' = {J u ^W°>. □ 

When w G d%), except in the case m = 1, only partial or local description of 5oo(w;) is available 
as shown in the next result and arguments that follow. However, as a simple corollary to Theorem 
13.41 5oo(u;) for w G 935 can be approximated by a limit of Soo(w n ) with {%} C 3°. To present 
the next result, we denote (v^ , . . . , Vi k ) by where = . . . , %k} C {1, . . . , m} and v G M m . 
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Lemma 5.9. For w G <9£>, there exists a nonempty index set M C {!,••■ ,m} such that the set 
{v(Q)m I f(^) e ^oo(^)} «s egua/ to the stable set of a system 

Vi = \vi l iex + ^2 AikVk + 9i( w )^ fori€M, 

keM 

which admits a unique equilibrium point. 

Proof. When w G d®, there exists some coordinate i such that rji(w) = Vi for any v G (£(u>). If not, 
then for each i = 1, . . . ,m there exists v G £(w) with Vi > i]i(w). A simple induction argument 
shows that A v + diagx(r](w)) has strictly negative diagonals in this case. This contradicts with the 
choice of w. The set M = {i G {1, • • • , m} \ rji(w) = Vi for any v G <B(w)} is, then, nonempty. 

We claim that Aij = for any i > j such that i £ A4 and j G .M c . Since is non- negative, it 
is enough to show that A^ cannot be strictly positive. Suppose A^ > for some i > j with i G M 
and j G A4 C . Since j G Ai c , there exists some u G £(w) with > rjk{w) for 7^ j and > r)j(w). 
Define 

In either case, it is easy to see V{ > rji(w) for i > j since > i]j(w). Also, the construction of 
gives Vi < Vi. This yields that the vector v := (v±, . . . ,Vi-i,Vi,Vi + ±, . . . ,v m ) lies in (£(u>) because 
f(v,w) < f(v,w) < 0. This is a contradiction to the assumption i G AL 

By construction, t){w)m is the unique equilibrium point of the subsystem iji = (l/2)y?Ij e x + 
Y,keM A ikVk + 9i(w), z G AL, of (]Ric-Vj) . For any y(0) G <Soo(w), z/ := lim t y(i) is an equilibrium 
point by Lemma [531 and it satisfies vm = t](w)m- Therefore, y(0)M belongs to the stable set of 
the aforementioned subsystem and the reverse inclusion is clear. □ 

We supplement Lemma 15 . 91 with an additional comment. When w G dD, every equilibrium point 
of (jRic-Vp is non-hyperbolic. In such cases, at least locally, the behavior of a solution trajectory is 
described by stable, unstable manifolds, and additionally, a center manif old, which is determined by 



a certain partial differential equation. We refer the reader to lPerkd 20011 ] for more deta ils about this 
topic. F or more detai 



2004J | or lOsinga et al. 



s on t heoretical and numerical studies of stable manifolds, see ICheng et al. 



20041 ] . Next, we prove the decomposition of S c 



Proof of Theorem \3.4\ If u = (v,w) G 5oo, then w G £>; otherwise, Corollay 15.71 implies that y(t) 
with y(0) = v blows up in finite time. Moreover, if u G S^, then we claim that v G S(w) with 
io G 33°. To see this, first note that ru G 5oo for some r > 1 sufficiently close to 1. Hence, rw G T> 
which is equivalent to <£(rw) 7^ 0. It then follows from Lemma 15.31 that <£°(w) 7^ 0. Therefore, 
w G £° and v G S^(w) = S(w). We then have C U w e»° 5 W x W- 

Conversely, pick any v G S(w) with w G 2)°. We want to show that, for any v! = (v',w ! ) 
sufficiently close to u = (v,w), we have u' G 5oo- Hence Uu;e£° ^( w ) x i w } — ^oo- ^° this end, we 
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consider, as in the proof of Lemma l5.8l a modified system r = {l/2)rj- + Ar + g(w') — g(w) where 
r(t) = y(t) — r](w), y(0) = v', and A = A v + diagx{r]{w)) . Recall that A has negative diagonals. 
This fact with straightforward calculations would imply that r(t) stays bounded above as long as 
\w' — w\ is small and r(0) is near the origin. Actually, it would be sufficient if we can find some to 
such that r(io) enters this neighborhood of the origin. However, linn, r(t) = when v' = v and the 
system for r(t) is smooth. Hence, r(i) continuously depends on its initial condition and thus such 
to can be found. 

Finally, noticing that dSoo = Uiogd Soo(w) x {w} we obtain the second statement from the 

first statement as well as from Lemmas 15.81 and 15.91 □ 



Now, we conclude this section with the proof of the characterization on exponential moments for 
affine processes. 

Proof of Corollary \3.6[ It follows from Theorem 12.11 that (|2 . 3[) is valid for any Xq as long as y(t) 
does not blow up by T. Then, the first statement is obtained. When E [exp (u 1 ^)] < oo for all 
T > 0, 



eMu T X T )] = \{uvy^u» + b D .u D + ±£b v - y(s) ds + Iy(T) • X% + V • X? . 

Lemma 15.51 implies that y(T) converges to an equilibrium point. Then (|3.2p follows from sending 
T — > oo in the above identity. □ 

6. Analysis of blow-up behavior 

In this section, we study solutions to (IRic-Vj) which blow up in finite time. We first introduce 
a change of variables in Section 16.11 Then we study the stability property of the system (|3.3|) and 
prove Lemma 13.81 in Section [6.21 Theorem 13.91 and Corollary 13.101 are proved at last. 

6.1. Blo w-up times. To s t udy t he blow-up time T*(u), we employ a change of variables investi- 
gated bv lElias and Gingoldl [20061 ]: 

(6.1) x(t) := 2y{t) 

l + yTT^WF 

where y(t) is a solution to (|Ric-V|) . This transform, which is equivalent to y = x/(l — \x\ 2 ) with 
\x\ < 1, maps W 71 onto the open unit ball in M. m . Moreover, \x(t)\ goes to the unit sphere whenever 
\y(t)\ goes to infinity. Therefore, this transform compactifies M m . Using this transform, we have 
the following representation of the blow-up time. 

Proposition 6.1. For each fixed u € W 71 x Ker A v , the blow-up time T*(u) is given by 

/•oo 

(6.2) T*(u) = / (l-R{s) 4 )ds, 

Jo 
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where R(s) := Y^iL\ x i{ s ) and x(s) solves 
(6.3) 



^ = (1 + R 2 )f(x, w) - 2{x ■ f(x, w))x, 



un 



th x(0) = 2v/(l + y/l + 4\v\ 2 ) and f(x, w) := (1 - R 2 ) 2 f (x/(l - R 2 ),w) . 



Proof. The proof of (|6.2p is essentially given in Section 3 of lElias and Gingoldl [2006J. We present 
their argument here for the reader's convenience. First, it is straightforward to check that (|6.ip 
yields 



(6.4) 



dx 
~dt 



1 



1-R(t) 2 



~ 2x-f(x,w) \ 



2v 



l + ^/l + 4\v\ 2 

Second, define s(i) via ds/dt = (1 - ^(t) 4 )" 1 for < t < T*(u) with R(t) = VEi^i^W 2 - Since 
R(t) 2 < 1 for t < T*(n), s(t) is a strictly increasing function with the unique inverse t(s). Let 
us write x(t(s)) and R(t(s)) as x(s) and R(s), respectively. Then, (|6.3|) follows from changing the 
variable t in (|6,4p to s. 
On the other hand, since 



d 



dx 



— (1 - = -2x • — = -2(x • f(x,w))(l - R( S ) 2 ), 

ds ds 

it follows that 1 — R(s) 2 = (1 — i?(0) 2 ) exp ^—2 J Q S x{r) ■ f(x(r),w)dr S j . Note that the integrand in 
this identity is uniformly bounded due to \x(r)\ < 1, then 1 — R(s) 2 > for any finite s. Thus, s 
maps [0,T*(u)) to [0, oo), hence (16. 2j) follows. □ 

The explosion time has the following property. 

Lemma 6.2. For each fixed u G W" 1 x Ker A v , define V{u) = {p G R + \ T*(pu) < oo}. ThenT*(-u) 
is a nonincreasing and differentiable function on V(u). 

Proof Note that ^ V because T*(0) = oo. Let us now prove the nonincreasing property of T*(-u). 
Without loss of generality, we assume 1 6 V{u). Then it suffices to prove T*(pu) > T*(u) for any 
p < 1. 

Consider a solution (y(t;p), z(t;p)) to (|2.2p with initial condition for p < 1. Then, y(t;p) 
satisfies = (l/2)y?Ij g i + (^4 v y)j + {p 2 /2)w T -KiW + p{A c w)i with y(0;p) = pv. Now, define 
y{t',p) = y(t;p)/p to obtain 



y = + A y + 2 



\ 



W T 7T m W 



( W T TTiW 



\ 



+ A c w, y{0;p) = v. 



+ A c w < X -yf + A v y + 
j \ w> ir m w J 

This implies that y(t;p) — f(y(t;p),w) < = y(t; 1) — f(y(t;l),w). It then follows from the 
comparison theorem before Lemma [A.4I that y(t;p) = py(t;p) < py(t; 1) for all < t < T*(u). This 
implies that T*(pu) > T*{u). 
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For the differentiability of T*(-u), first observe 



-f^(s) 2 = , 4 — — x [^yMs)-p)^- yi {t{s)- P ) 

dp V 1 + My(t(s); P )\ 2 (l + Vi + My(t(s); P )\^ v - dp 



a=l 

where t(s) is the function defined in Proposition 16.11 The analysis in the previous paragraph 
implies that (d/dp)(y(t;p)/p) > 0, which in turn gives dy/dp > y/p. On the other hand, Remark 1 5. 6 1 
implies that lim t ^. T *^ pu ^ yi(t',p) = oo. Here i is one component such that yi(t;p) explodes at T*(pu). 
Combining the previous observations, we obtain that 

yi{t{s);p)^j-yi{t{s);p) > ^yf(t(s);p), 

up p 

which is unbounded from above. Moreover, thanks to Lemma l5.5l yi{d/dp) yi dominates other non- 
explosive components. Therefore, there exists some sufficient large sq such that (d/dp)R(s) 2 > 
for all s > sq. Hence it follows from (|6.2h that 

d d / f s ° f°° \ r°° d 

^-T*(pu) = -£-(/ (1 - R(s) A ) ds+ (1 - R(s) 4 ) ds)=- 2R(s) 2 j-R(s) 2 ds, 
dp dp \J J sq J J dp 

where Fubini's and Tonelli's theorems are applied to integrals on [0, sq] and [sq, oo), respectively. □ 



Not only does the change of variable (16. If) provide an expression for blow-up times, but also it 
helps to study the blow-up rate of solutions. 

Lemma 6.3. Suppose that T*(u) < oo for some u G W 71 x Ker AP , and that I < m is the first 
component of y(t) that blows up at T*(u). Then, ]hn.^*r u \(T*(u) — t)yi(t) = c for some positive 
constant c. 

Proof. For notational convenience, we write T* for T*(u). First, we observe that I must be in I. 
Otherwise, y t = A u yi + Y,k=i A ikVk + 9i(w) where yi(t), . . . , y t -i(t) are finite on [0, T*], from which 
we infer that yi(t) must be finite in [0, T*]. This contradicts with the choice of I. 

Now, we apply the compactification (I6.ip to (yi, . . . ,yi). The resulting function x(s), which is a 
vector valued function of length I, satisfies (16. 3ft . It then follows from the choice of I and (16. If) that 

, , v® f° i<l . 



The vector e; can be easily verified to be an equilibrium point of f|6.3f) . If we denote the right hand 
side of (|6,3p by h(x) = (hi, . . . ,hi), then 

^ = 2 Xj h + (1 + R 2 )§^ ~ 2 *kfk) Sii - *H U + E *kj£) , 



where 5ij is the Kronecker delta and 



^ = Xiligi^- - 2xj E + (1 - &)A*j ~ 4 (! " R 2 )x j g(w). 



k=l 



ON EXPONENTIAL MOMENTS OF AFFINE DIFFUSIONS 29 

From above calculations, we obtain that the Jacobian matrix of h at e/ is —I\. It is clear that 
the eigenvalues Ai, . . . , A; of this Jacobian matrix are non-resonant, i.e., there is no (mi, . . . ,mi) 
such that m k £ {0| U N , Ylk=i m k — 2, and A& = J2\=i m -i^i f° r k G {!>•••) 0- Theorem 4.1 in 



Elias and Gingoldl 20061 ] now implies that y^ =1 Vk(t) 2 ~ c(T* — t) 1 as i f ^* f° r some positive 
constant c. Since is the first component in y(t) that explodes at T*, we consequently have 
yi (t) ~ c(T* - i)- 1 , as t T T*. □ 

6.2. Blow-up regions. Before proving Theorem 13.91 and Corollary 13.101 in this subsection, let us 
present a stability property of the system (|3.3p and prove Lemma [ 



Lemma 6.4. If the trajectory of 13. 3\) is bounded, then it converges to an equilibrium point. If the 
trajectory is unbounded, then it blows up in finite time. 

Proof. Since x m+ \ clearly converges to zero, which is the (m + l)-th coordinate of any equilibrium 
point of (|3.3p . we only need to prove the first statement for the first m coordinates. To this 
end, we prove by an induction on i. But we present the induction step only. The case i = 1 is 
straightforward. 

Now suppose that x±, . . . , converge to the first i — 1 components of some equilibrium point. 
Recall that x { satisfies x { = ^xfl^x + {TAue~ s - 1) Xi + T J2l~=i Ake~ s x k + Te~ 2s g(w). If i £1, 
it then follows from the induction assumption and Lemma lA.21 that lim^oo Xi(s) = 0. If i G I, 
define x(s) := Txi(s) + TAae~ s — 1, then x satisfies 



i-l 

= 2** + T2 Y1 A ike~ s x k + T 2 e- 2s g(w) - (TA^ 5 - l) 2 - TA u e 



1 z2 , t-,2 a „-s_ , t^2„-2s„/„..\ 1 (, 

k=l 



Lemma I A . 1 1 implies that i(oo) := liiaig^.^ x(s) E M and (l/2)x(oo) 2 — 1/2 = 0. Therefore, x(oo) = 
±1, which implies Xj(oo) = or 2/T. Hence, in the above two cases, Xj(s) converges to the i-th 
coordinate of some equilibrium point. This concludes the induction step. 

For the second statement, let us assume that Xi(s) is the first component whose trajectory is 
unbounded. It is necessarily that i G I; otherwise, we can utilize Lemma IA.2I to deduce that 
Xi(s) — > 0. Moreover, the trajectory of Xi(s) is bounded from below. Indeed, using the comparison 
theorem for scalar ODEs, we can see that Xi(s) > x(s) where x solves 

x = (TAue' 3 - l)x + T^A lk e- s x k {s) + T e - 2s gi (w), x(0) = Xi (0). 

l 

Again Lemma lA.21 implies that x is bounded. Hence, the trajectory of Xj is bounded from below. 
Now, since x\, . . . , are bounded by assumption, we can find a positive constant C such that 

T 

Xi > —x 2 + (TAae~ s - l)x; - C =: h(xi, s). 

Moreover, there exists a sufficiently large r such that r > T _1 (l + \/l + 2TC) and h{r, s) > for 
any s > 0. On the other hand, the trajectory of X{ is unbounded from above by assumption, but 
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bounded from below. Therefore, we can find some so such that Xj(so) = r, and then Xi(s) is strictly 
increasing from so onward. As a result, 

T 2 

Xi > —x { —Xi — C, S> Sq. 

However, notice that Xi(so) > T _1 (l + yl + 2TC) which is the unstable equilibrium of the previous 
ODE. Now it is immediate to check that Xi(s) blows up in finite time. The second statement is 
proved. □ 

Proof of Lemma \3.8[ To show the continuity of T*(-) on V, let us denote the solution to (IRic-VI) by 
y(t;u) where its dependence on u is explicitly indicated. Let T*(u) = T. The associated solution 
x(s;u) to (|3.3p does not blow up in finite time; otherwise, y(t;u) explodes before T. Lemma 16.41 
implies that the trajectory of x(s;u) is bounded, which in turn yields that lim s _ s>00 x(s; u) = v' for 
some equilibrium point v' of (|3.3p and u' = (u, 0) as set in the subsequent discussion. 

However, this point has at least one nonzero coordinate. This is because, for the first index I 
such that yi(t; u) blows up at T, we have from Lemma 16.31 that lim i ^(T — t)yi(t; u) = c > 0, which 
implies 

(6.5) lim xi(s;u) = \imxi ( — log ( 1 — — J ;u J = lim ( 1 — — J yi(t;u) = — . 



t\T \ °\ TJ J trr\ T ) y ' T 

Combined with the characterization of equilibria for (|3.3p . it is necessarily that c = 2; see the 
discussion after (|3.3|) . Moreover, Uj = if j ^ I. From this blow-up behavior of y(t;u), it is 
easily deduced that there exists some i £ I such that yi/\y\ converges to a positive constant while 
lirnt yj/\y\ = for all j g 1. 

Now, let us consider z(s;u) a solution to (|6.3p with z(0;u) = 2v/(l + yl + 4|v| 2 ). From zj\z\ = 
y/\y\ and lim s \z(s; u)\ = 1, lim s 2j(s;n) > and lim s Zj(s; u) = for all j ^ X. Therefore, for 
some positive C, small 5 and sufficiently large sq, ^ fcG j Zk(s;u) 2 > C and Zi(s;u) > — 5 for all 



Lefschetz 



I = 1, . . . , m and s > sq. Furthermore, since z(s; u) continuously depends on u (see e.g. 
19631 ] ). there exists a sufficiently small neighborhood U of u, such that ^fcez z k( s ', u ') 2 > C and 



zi(s; u') > —5 for alH = 1, . . . , m, s > so, and all u' G U. 

Thanks to the analysis in the last paragraph, we can find a sufficiently small e > such that 
Sfcex z k( s 'i u ') > 6 f° r an s — s o anc ^ u> € U. As a consequence, for an even larger so, we can see 
that 

z(s;v!) ■ f(z(s;u'),w) 

= \Y. z ^ u 'f R ^ u ') 2 ) z ^ u) T A v z(s; v!) + (1 - R(s; u'f) 2 z(s; u') T g(w) > | 

fcez 
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for all s > so an d u' £ U because the second and third terms become small as \z\ < 1, and R 2 
converges to 1 as s increases. Therefore, eventually we get 

1 - R{s; u'f < (1 - R(0) 2 ) exp ^-2 J ° z(r; u) ■ f(z{r; u'),w)dr - e(s - s ] 

using the functional form of R(s; u) 2 in the proof of Proposition 16. 11 This facilitates the application 
of the dominated convergence theorem to conclude that lim n T*(u n ) = J °° lim n (l — i?(s;u n ) 4 ) ds 
for a sequence of initial conditions u n E U that converges to u. The right hand side of the previous 
identity, then, is equal to T*(u) due to the continuous dependence of x(s;u) on u. □ 

Proof of Theorem \3.9l Let us consider the case of OSt first. The beginning paragraph of the proof of 
Lemma argues that the solution y(t) of (|Ric-Vp that explodes at T is associated with the function 
x(s), the solution to (|3.3p with x(0) = (y(0), 1), and the limit lim s x(s) is equal to some nonzero 
equilibrium point v' = (u, 0) of (|3.3p . In other words, y(0) £ W£(w, T). Hence, {v \ T*(v, w) = T} C 
U^ WZ( W ,T). 

For the converse, suppose v S W£(w,T) for some nonzero equilibrium point v' = (^,0). Recall 
the discussion following (|3.3p . and thus Vi is either zero or 2/T. Since lims-j.oo x(s) = (^,0) when 
x(s) is the solution to (|3.3p with x(0) = (v, 1), the same computation for Xi(s) as in (|6.5p yields 
that y(t) is finite for all t < T and lim^T^T — t)yi{t) = 2 for any i such that Vi ^ 0. Therefore, 
v € {v | T*(v,w) = T}. This completes the proof of the first assertion. The second assertion then 
clearly follows. 

As for it is already noted in Section [3.21 that \im s ^ 00 x(s) = when y(-) blows up after T. 
Hence, C U^eRer^ 17 ^ / o( s '^) x i w } ^ s c l ear - The reverse inclusion also easily follows: If x(s) 
converges to zero, then T*(u) > T, but T*(u) = T cannot happen because, otherwise, the limit of 
x(-) would be a nonzero equilibrium point as shown above. □ 

Proof of Corollary \3.10l The first statement follows from the similar reasoning as in Corollary 13.61 
Hence, we focus on (|3.4p . From the transform formula (|2.3p . for S < T where u G dSx, 

logE exp (u T X s ^j = |(n D ) T vr n D + S b v ■ u v + b v ■ y(s)ds + y(S) ■ X% + u v ■ Xg. 

Therefore, lim stT (T-5) logE[exp(n T X s )] = lim sr r(T-S) f Q S b v ■y(s)ds + lim stT (T - S)y(S) ■ X^ . 
But, we know from Theorem 13.91 and (|6.5p that lim,g^T(T — S)y(S) = Tv for some z/, the first m 
components of some nonzero equilibrium point of (|3,3p . Then, it is a simple exercise to show that 
(T — S) Jq b v ■ y(s)ds converges to zero. Now, the result is immediate. □ 

7. Conclusion 

In this paper, we study the long-term and blow-up behaviors of E [exp (u t Xt)] for multi- 
dimensional affine diffusion X with some hierarchical structure between components. Analyzing 
solution behaviors of a multi-dimensional Riccati system, which is associated with a given affine 
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diffusion process via the transform formula, we completely characterize sets of u such that exponen- 
tial moments are finite for all time or only before a fixed time. These sets are decomposed into the 
unions of stable sets for equilibrium points of the Riccati system or its transformed version. Then, 
we compute certain limits of exponential moments which provide detailed descriptions of behaviors 
of affine diffusions. 

When the log-return of discounted stock prices is modeled by a linear transformation of affine dif- 
fusion processes, our results identify the long-term and blow-up behaviors of stock prices, especially 
in the case where the stock price moment is not explicitly known. These results provide a handle 
to investigate the implied volatility asymptotics for large-time-to-maturity, deep-out-of-money and 
deep-in-the-money options. We presented several examples to illustrate this point. Theoretically 
and practically, it still remains an interesting topic to extend the analysis of this paper to affine 
processes with jumps or affine processes on more general state spaces. 

As a final remark, it is well known that in the literature of affine processes bond options and some 
other fixed income products can also be expressed in semi-closed form using the Fourier inversion 
formula. As long as the long-term growth rate of the underlying process satisfies the assumptions 
of the Gartner-Ellis Theorem, we can calculate the asymptotic behaviors of the price of such a 
product, which are possibly useful in obtaining quantities that are analogues of the Black-Scholes 
implied volatility. We leave this as a potential future research topic. 

Appendix A. Auxiliary results on ODEs and matrices 

Lemma A.l. Let us consider a scalar ODE x(t) = x{t) 2 /2 + h(t) with h(oo) := lim^oo h(t) 6 R. 
If the whole trajectory {x(t) : t > 0} is bounded, then h(po) < 0, x(oo) := lim^oo xit) 6 R, and 
x(oo) 2 /2 + h(oo) = 0. 

Proof. Let us prove h(oo) < first. Otherwise, h(oo) > 25 for some positive constant 5. As a 
result, there exists some to such that x(t) 2 /2 + hit) > x(t) 2 /2 + 5, for a l l t > tp. It follows from the 



comparison theorem for scalar ODEs (see e.g. Chapter II of iHartmanl [19821 ]) that x(t) > y(t) for 
any t > to, where y(t) is a solution to y = y 2 /2 + 5 with y(to) = x(to)- However, a simple analysis 
of the previous ODE yields that yit) blows up in finite time. This contradicts to the assumption 
that the trajectory of x(t) is bounded. 



To prove the rest of the statements, we shall first show that limsup^oo xit) < \J — 2/i(oo). If 
not, then there exists 5 > such that limsup t ^ 00 x(t) > \J '— 2/i(oo) + 25. Then, we can find to such 
that x(fo) > \J — 2/i(oo) + 25 and h{t) > h(oo) — 5 for all t > to- Next consider y_(t), which satisfies 
y_ = y 2 _ /2 + /i(oo) — 5 with y~(to) = x(to). The comparison theorem implies that x(t) > y~(t) for 
t > to- However, y_ explodes to infinity because y(to) > \/—2h(oo) + 25 and because the value on 
the right hand side is the unstable equilibrium point of the ODE satisfied by y_. This contradicts 
to the boundedness assumption on {x(t) : t > 0}. 
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Now, if lim inf t-^oo x(t) > y— 2/1(00), then, combined with the result from the last paragraph, 
we have lim^oo x(t) = \f—2h(oo) and we are done. To deal with the other case liminf^oo x(t) < 
■sj — 2/i(oo), we separate h(oo) = and h(oo) < cases. 

When h(oo) = 0, if lim int^oo x(t) < y / '— 2/i(oo), then we have lim inf^oo x(t) < V25 for any 5 > 
0. Thus, there exists to such that x(to) < \^25 and h(t) > —5 for all t > to- Hence, lim^oo y_(t) = 
—y/25 because y_(to) = ^(^o) is less than the unstable equilibrium of y_. Combining this with x(t) > 
y~(t) for t > to, we obtain liminf^oo x(t) > linit-j-oo U-(t) = — v25, which implies liminf^oo x(t) > 
thanks to the arbitrary choice of 5. This is a contradiction. 



When h(oo) < 0, if liminfi„ >00 x(t) < a/— 2/i(oo), then there exists a sufficiently small positive 5, 
such that /i(oo) + 5 < and liminf^oo x(t) < '— 2/i(oo) — 2<5. We can find to such that x(to) < 
y 7 — 2/i(oo) — 25 and /i(oo) — 5 < /i(t) < h(oo) + 5 for all t > to- Consider y + which satisfies y+ = 
y\/2 + h(oo) + 5 and y+(to) = ^(to)- Note that ^/—2h{co) — 25 is the unstable equilibrium of y + , we 
then have from the comparison theorem that — ^—2/1(00) + 25 = limt_ s . 0O y~(t) < liminft_ i , 00 x(t) < 
limsup^oo x(t) < lim^oo y+(t) = — y/—2h(oo) — 25. Since 5 can be made arbitrarily small, we 
conclude from previous inequalities that lim^oo x(t) = — y— 2/i(oo). □ 

Lemma A. 2. Let us consider a scalar ODE x(t) = (ae~ l — l)x(t) + g{t) with g € C 1 and g(oo) := 
limt_>.oo d(t) S R. Then, x(oo) := linit-s-oo x(t) G R and x(oo) = 5(00). 

Proof. Consider a new function y(t) = exp (ae _t + t) x(t). Then, y satisfies y(t) = exp (ae~* + t) <?(t) 
and consequently, 

x(t) = exp (a — ae _t — t) x(0) + exp (— ae~ l — t) / exp (ae _s + s) g(s)ds. 

Jo 

Choose an arbitrary e > 0. Then, we can find a large T = T(e) > such that \g(t) — g(oo) \ < e 



for all t > T and 



1 - e ae T 



< e. Next, we compute for t > T 



f e ae " +s g(s)ds 







e - * / e ae s+s g{s)ds + e- t I e ae s+s g(oo)ds + e~* / e ae s+s ( 5 (s) - 5(00)) ds. 
io Vt Jt 



The last term is bounded by 



e-* l t e ae -' +s (g(s)-g(^))ds 



T 



-T 



1 (l- e ^)< e (l + e )(l-e^) 



using < e s < e T . For the second term, we obtain e ae T (l — e T ') < e * J^e ae " +s ds < 
(l — e T_t ) if a < (inequalities are reversed if a > 0). Therefore, by the assumption on T, 



5(00)6"* / e ae ~ s+s ds - g(oo) (1 - e T ~ l ) 



<e| 5 (oo)|(l-e T - f ) 



These calculations yield lim^oo e 4 J * e ae " +s g{s)ds = g(oo) + c with |c| < e (|5(oo)| + 1 + e). Since 
e is arbitrary, we can conclude that lim^oo x{t) = 5(00). □ 
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The following definition and property can be found in Chapter 6 of lBerman and Plemmons 



Definition A. 3. A square matrix is called a nonsingular M-matrix if it has non-positive off- 
diagonals and every real eigenvalue is positive. 

If M is a nonsingular M-matrix, then M^ 1 > for all i,j. Before we present the next result, we 
recall the following extension on the comparison theorem for scalar ODEs. A function / : R m — > R m 
is quasi-monotone increasing, if fk{x) < fk(y) for any x, y such that x k = y k for some k and Xj < y-j 
for any j ^ k. Suppose that / is quasi-monotone increasing and locally Lipschitz, then for any two 
differentiable functions x,y : R+ — > M m with x(a) < y(a), 

x(t)-f(x(t))<y(t)-f(y(t)), Vt>0 => x{t)<y{t), Vt > 0. 



A proof can be found in IVolkmannl [1972 ] 



Lemma A. 4. Let 5 € (0, 1) and M be a nonsingular M-matrix. Consider two solutions x(t) and 

r<y\ 

y(t) of x = (l/2)xj — Mx with x(0) = 5y(0) and I C {l,...,m}. Then x(t) < 5y(t) whenever 
they exist. 

(2) 

Proof. We define / via f(x) := (l/2)xj- — Mx. It is clearly locally Lipschitz as well as quasi- 
monotone increasing. Indeed, for any x < y with xt = yk for some k and Xj < yj for all j ^ k, we 
have f k (x) = (l/2)xp fce i - XX=i M kjXj < (l/2)xp fee i - M kk x k - J2j^k M kjVj = fk(y), thanks to 
M kj < for j / k. 

Now consider z(t) := x(t)/5 which satisfies i = (5/2)z^ — Mz with z(0) = y(0). Then, we have 

i(t) - f{z{t)) = ^4 2) < = y(t) - f(y(t)). 
We conclude z(t) < y{t) from the above comparison theorem. □ 
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